Geometrical Pump in Quantum Transport: Quantum Master Equation Approach 
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For an open quantum system, we investigate the pump current that is induced by a slow modula- 
tion of control parameters on the bases of the quantum master equation and full counting statistics. 
We find that the average and the cumulant generating function of the pump are characterized by the 
geometrical Berry-phase-like quantities in the parameter space, which is associated with the gener- 
ator of the master equation. By our formulation, we can discuss the geometrical pump under the 
control of the chemical potentials and temperatures of reservoirs. We demonstrate the formulation 
by spinless electrons in coupled quantum dots. We show that the geometrical pump is prohibited 
for the case of non-interacting electrons if we modulate only temperatures and chemical potentials 
of reservoirs, while the geometrical pump occurs in the presence of an interaction between electrons. 

PACS numbers: 05.60.Gg, 72.10.Bg, 73.23.-b, 73.63.Kv 
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I. INTRODUCTION 

When a quantum system is slowly and periodically 
modulated by two or more control parameters such as 
gate voltages, a net number of particles can be trans- 
ported per period of the modulation even in the ab- 
sence of dc driving force (e.g., bias voltage). This phe- 
nomenon is known as an adiabatic pump. The original 
idea of the adiabatic pump was proposed by Thouless 1 , 
where a pump current of a closed quantum system is 
related to the Berry phased of the ground state of the 
Hamiltonianii 3 -! 4 -. Since then, experimental and theoreti- 
cal studies of the adiabatic pump have been developed 5-7 . 

Later, the idea of the adiabatic pump has been ap- 
plied to open quantum systems both in experiments 8 
and theories^—. Theoretically, the average pump cur- 
rent can be expressed by the Berry phase associated 
with the scattering matrisi^. The cumulant generat- 
ing function of pump current can also be described by 
geometrical quantities^ 2 -. For this reason, the adiabatic 
pump is referred to as the geometrical pump. The above- 
mentioned theories are based on the time-dependent scat- 
tering theory^ 3 -. 

Similar phenomena have been studied in stochastic sys- 
tems described by classical master equationls^I, which is 
referred to as the adiabatic stochastic pump. The pump 
current in the classical stochastic pump has also geomet- 
rical properties; the cumulant generating function of the 
pump current is expressed by a Berry-phase-like quan- 
tity that is associated with the generator of the classical 
master equation. We shall refer to this quantity as the 
Berry-Sinitsyn-Nemenman (BSN) phas o 16 ' 17 . There have 
also been several works on the adiabatic pump for quan- 
tum open systems described by quantum master equa- 
tion (QME) 2 ^— . However, the geometrical effects in the 
adiabatic pump described by QME have not been fully 
understood. 

In this paper, we investigate the quantum adiabatic 



pump on the basis of QME, and derive general formulae 
of the cumulant generating function and average of the 
pump. These formulae are geometrical and expressed by 
a quantum analogous of the BSN phase which is associ- 
ated with the generator of the QME. Since our theoretical 
framework is independent of the details of the system, it 
is useful for analyzing a large variety of applications of 
the adiabatic pump such as a qubit rotation in quan- 
tum dots and a spin pumping. We note that the master 
equation approach is suitable to analyze systems that in- 
clude interaction between particles^ 4 -, dissipations, and 
decoherences 2 ^. 

This paper is organized as follows. In Sec. QI] A, the 
QME approach for an open system is discussed. The full 
counting statistics in the QME approach is discussed in 
Sec. HII B. For the adiabatic pump, the geometrical formu- 
lae (BSN phase expressions) for the cumulant generating 
function and the average of the pumped quantity are de- 
rived in Sec. |TT]C. In this formulation, the temperatures 
and chemical potentials of reservoirs are the control pa- 
rameters. This is in contrast to most of the conventional 
studies on the adiabatic pump, where only the parame- 
ters in the system Hamiltonian or in the coupling with 
the reservoirs are considered. In Sec. we demonstrate 
our theory by the spinless electron transport in quantum 
dot systems. For non-interacting cases with (Sec. IIIII A) 
and without (Sec. lIIll B) the rotating wave approximation 
(RWA), we find that no geometrical pump occurs if we 
control only the temperatures and chemical potentials of 
reservoirs. In contrast, the geometrical pump occurs for 
an interacting electron system even when we control only 
the reservoir parameters (Sec. IIIII C). Section ITVl is de- 
voted to the discussion and conclusion. In Appendix [XJ 
we describe the details of the QME in the framework of 
the full counting statistics. In Appendix [B] we show the 
detailed derivation of Eq. (flTf . In Appendix [Cj we ver- 
ify the equivalence between without and within the RWA 
concerning the unit-time generating function of the cur- 
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Reservoir L System S Reservoir R 

FIG. 1. (Color online) Illustration of our setup with two reser- 
voirs (L and R). 



rent between the system and reservoirs in a steady state. 

II. GENERAL RESULTS 
A. Setup 



where A„(a) is an eigenvalue of /C, n is a label of the 
eigenvalues, and p n (p) is the corresponding right eigen- 
vector. 

By introducing the Hilbert-Schmidt inner product of 
linear operators A and B on TLs as Trs(ATS), we de- 
fine the adjoint lO of the QME generator such that 
Tr s [(£ti)tf?] = Tr s (lt/CB) holds for any A, B. We 
then have the left eigenvalue equation for /C: 

K)l n {oc) = \* n (<x)L(a), (4) 

where £ n (a) is the left eigenvector corresponding to the 
eigenvalue X n (a). In the following, we assume that K, 
has the zero eigenvalue Ao = without degeneracy, so 
that ICpa = and K,Ho — hold. This implies that the 
QME has a unique steady solution p ss = po for fixed ex. 
We note that £o(ot) — 1 (identity operator) holds for any 
a. 



We consider a quantum system S that is weakly cou- 
pled to reservoirs {R&}, where 6 is an index of reservoirs 
(see Fig.[T]for a schematic). The total Hamiltonian of the 
coupled system is H tot = H s + J2b( H b + H S b), where H s 
is the system Hamiltonian, H b is the Hamiltonian of the 
6th reservoir R;,, and H$ b is the interaction Hamiltonian 
between S and R^. If the interaction between the system 
S and the reservoirs is weak, the dynamics of S can be 
described by a QME for the reduced density matrix of 
S, which is denoted as p. Suppose that the initial state 
of the system S is decoupled with the reservoirs. Then, 
up to the second order in the system-reservoir coupling 
(Born approximation) with the Markov approximation 26 , 
the QME for the system S reads 



dp{t) 
dt 



JCp(t), 



(1) 



where 



Kp=^[H s ,p} + ^2v b p, (2) 

b 

1 f°° 

V bP =-— J dt'T n [H sb ,[H sb (-t'),p®p b }]. 

Here the symbol " " " stands for the interaction picture 
with respect to Hs + ~^Z b H b , Trb represents the trace 



over the 6th reservoir, and p b 



,-Pb{H h -n b N h ) 



jZ b is the 



grandcanonical distribution with the inverse temperature 
j3 b , chemical potential p b , and the particle number opera- 
tor N b of the 6th reservoir. The time-evolution generator 
JC of the QME depends on several parameters; the sys- 
tem parameters in Hg and Hs b such as the energy levels 
of quantum dots and the tunnel barriers between them, 
and the reservoir parameters, {/?&} and {p b }. We write 
the set of these parameters as a. The right eigenvalue 
equation for K, is written as 



B. Full counting statistics 

We consider the statistics of a quantity Aq transferred 
from the reservoirs to the system S during a time interval 
t. The measurement scheme of Aq is as follows. First, at 
t = we perform a projection measurement of a reservoir 
variable Q to obtain a measurement outcome go- We as- 
sume [Q, N b ] = and [Q, H b ] = for any 6. For t > 0, the 
system S undergoes the time evolution with interacting 
with the reservoirs. At t = t we again perform a projec- 
tion measurement of Q to obtain another measurement 
outcome q T . Then Aq is defined as Aq — q T — qo . 

The cumulant generating function of the statistics is 
given by S T (x) = In / P T {Aq)e l * A idAq, where P T {Aq) is 
the probability of Aq during r. \ is called the counting 
field, and the derivatives of S T (x) give the cumulants of 
Pr{Aq); e.g., (Aq) T = dS T ( X )/d(i X )\ x =o. Note that if 
Q is the 6th reservoir's particle number N b (Hamiltonian 
H b ), then {Aq) T /r is the average of the particle (energy) 
current from the 6th reservoir into the system S. 

For calculating S T (x), we employ a method devel- 
oped in the context of the full counting statistics 2 ^. In 
this method <SV(x) i s obtained from the solution of the 
modified equation of motion which is governed by the 
X-depcndcnt Hamiltonian. In the QME approach, this 
reads <SV(x) — hiTrs/5 x (r), where Trs is the trace over 
the system S and p x (r) is a solution of the generalized 
QME (GQME): 



dp x (t) 



= W- 



dt " xh 

Here the modified generator IC X is given by 



bP ' 



(5) 



(6) 



Kp n (a) = \ n (a)p n (a), 



(3) 



h 2 



dt'Tr b [H Sb ,[Hsb(-t'),p®Pbh 
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P0- x and O x 



where [0,P] X = 6 X P 
e -ixQ/2() e ixQ/2. See Appendix |A] for details. 

We denote the left and right eigenvectors of K, x 
(for fixed a) corresponding to the eigenvalue A*(cti) as 
£ x (a) and p*(a), respectively. They are normalized as 
Trs(^mPn) = $mn- We assign the label for the eigenvalue 
with maximum real part to n = 0. Then p x (t) ~ e A ° r 
for large r, which results in lim T _). 00 SV(x)/ r = Ag. Thus 
Aq(cx) is the unit-time cumulant generating function of 
the steady state for fixed a£± Note that if we set x = 0, 
the GQME reduces to the original QME, and £ x and p x 
to 1 and pss, respectively. 



C. Geometrical pump 

We slowly modulate the parameters a along a curve C 
in the parameter space during a time interval r. If the 
system is in the instantaneous steady state for the value 
of ctt at each time t in the whole of the process, the cumu- 
lant generating function for Aq for this process is equal 
to the time integral of the unit-time cumulant generat- 
ing function Aq (ctt) of the instantaneous steady state. In 
general, however, there exists additional (pumped) con- 
tribution; S T {x) = Io dtX o( a t) + St*{x)- We call the 
latter contribution the excess part. The excess part is 
intrinsic in the transitions between the steady states, and 
is of our interest. 

We now derive the geometrical expression of the excess 
part of the generating function by using the method sim- 
ilar to those in Refs^ 6 -^. First, to solve the GQME for 
a given curve C of a in the parameter space, we expand 
p x (t) as 



P x (t) - 



c n (t)e A ^p x (a t ) 



(7) 



where A%(t) = f Q dt' A*(ctt<) . Inserting this equation into 
Eq. ((SI) and taking the Hilbert-Schmidt inner product 
with l x (a t ), we obtain 



dcpjt) 
dt 



5>(t)e A ^Tr s ^(a 



dt 



(8) 



We assume that the modulation of the parameters a is 
sufficiently slow. Then we can approximate the sum on 
the right hand side of Eq. ((5J by the contribution only 
from the term with n = 0. By solving this adiabatic 
approximation equation we obtain 



co(t) = c (0)exp 



- / Tr s (i x Ha)dp x (a)) 



(9) 



where dp* {a) = da. ■ dp x (a) / da. If the initial state 
of the system S is p x (0) = p S s(&o), then Co(0) — 

Tr s in (a )p ss (a ) . We again use the adiabatic ap- 



proximation to obtain 

p x (r) ~co(r)e A oM^(a T ) 

= e A o^)^(a T )Tr s (^ t (a ) / 5ss(ao) 

x exp 



$t(a)dp*(a)) 



(10) 



Thus we obtain the excess cumulant generating function 
S° K (x) = S T (x) — Aq(t) for the slow modulation: 



S?(x) = - [ Tr s (l x \a)dp x (a) 



+ In Tr s 



(^(oojA-Coo)) +lnTr s ^(ct T ). ( 



11) 



The right-hand side of Eq. (fTTj) is analogous to the Berry 
phase, which depends only on the line integral along the 
curve C (not on time). We note that Aq(t) corresponds 
to the dynamical phase. 

By differentiating Eq. (|11[) with respect to ix, we ob- 
tain a geometrical expression of the average excess in the 
quantity Aq: 



Tr s ( #(a)^(a) )■ ,/«. (12) 



where i' = dl x / d{ix)\ x =o- Equations <[TTJ) and (fl2|) 
are regarded as quantum versions of the Berry-Sinitsyn- 
Nemenman (BSN) phases for the cumulant generating 
function and the average, respectively, in slow paramet- 
ric modulatio n 16 ! 21 . We refer to Trs (£q 9/5 ss /9a) as the 
BSN vector potential for average excess quantity. 

Equality (|12l) implies that a finite net quantity of Aq 
can be transferred to the system S for a slow cyclic modu- 
lation of the parameters (i.e., for the case where the curve 
C is a closed loop) even if there is no dc driving force such 
as temperature and chemical potential differences. This 
is an adiabatic pump. By the Stokes theorem, Eq. (fT2|) 
is rewritten as 



(Aq)™ 



dam A du r . 



(13) 



for a cyclic process, where A is the wedge product, Sq is 
a surface enclosed by C, and 



Tr s 



dip dp s; 
dam da, 



dip dp s 
da n da r , 



(14) 



We refer to Eq. (|14l) as the BSN curvature for average 
pumped quantity. 

For some cyclic processes, the excess pump (Aq)° x van- 
ishes. A sufficient condition for the "no-pumping" is that 
Fa m a n = holds in Sc for all (a m ,a n ). We note that 
if the whole of a curve C (not necessarily closed) lies in 
a region of the no-pumping condition, then the average 
excess pump does not depend on the whole of C but de- 
pends on only the initial and final points of C. 
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III. APPLICATION TO SPINLESS ELECTRON 
TRANSPORT IN QUANTUM DOT 

In this section, we apply the general framework ob- 
tained in the previous section to the transport of spinless 
electrons in quantum dot systems connected to two elec- 
tron reservoirs (6 = L, R). The spinless electron model 
can be used for analzying phenomena in which the spin 
degrees of freedom are irrelevant. The reservoir Hamil- 
tonian is given by Hb = J2k ^bkc\ k Cbk- Here hilbk is 
the fcth mode energy of the electron in the 6th reser- 
voir, and &l k (dbk) is the corresponding creation (anni- 
hilation) operator, which satisfies {cl k ,c b 'k'} — <$fcfc'<W, 

{^bk'^l'k'} = {cbk,Cb'k'} = 0. We here take the quantity 
to be counted as the electron number transferred from the 
reservoir L to the system S, i.e., Q = Nl = J2k ^LkpLk- 



A. Non-interacting electron model with RWA 

First we consider a model of non-interacting spinless 
electrons in a series of N single-level quantum dots cou- 
pled to two reservoirs (6 = L,R). The Hamiltonian of 
the system S in this model is 

H s = £ ^U* + 12 + h - c -' ( 15 ) 

where £j is the level energy of the ith dot, d] (di) is the 
creation (annihilation) operator of the electron in the zth 
dot, and v is the transfer probability amplitude between 
the dots. This Hamiltonian can be diagonalized by a 
unitary transformation as H§ = y\ tuj-ja^aj. Here fkjjj 
is the jth mode energy of the electron in the system S, 
and at is the corresponding creation (annihilation) 
operator of the electron, which satisfies the canonical 
anti-commutation relations of fermion: {at, a.,'} = Sjf, 

{at, at,} = {oj,Ojv} = 0. The coupling Hamiltonian 
between the system S and the 6th reservoir is given by 
Hsb = J2j,k V bjka]c bk + h.c. 

To derive the GQME for this system we use the 
Born (up to the second order in Hgb) and the Markov 
approximations^ 6 -. Furthermore, we adopt the RWA, 
which is a coarse-graining of the time evolution on the 
time scale longer than that of the system evolution with- 
out the coupling with the reservoir s 26 ' 27 . For x = the 
RWA leads a Lindblad form of the QME and guarantees 
the complete positivity of the time evolution. We carry 
out the RWA by averaging over the rapidly oscillating 
terms in the Born-Markov GQME in the interaction pic- 
ture (see Appendix lAl for the detail). Then we obtain the 
generator /C x of the GQME in the form of JC X — J2j ^x,j- 
Here the GQME generator fC x j for the jth mode is given 



by22 

-tffi 12 (% J {u j ){a)a ] p + pa)a ~2e~^afpa]} 

b=L,R 

+ ^(wiHojfitp + pa 3 a) - 2e*»ot / 3aj}) , (16) 

where H^ h ee - *+(w i )a i a]}/2fi 

is the Lamb shift Hamiltonian, $^(w) = 
27r^ fe \Vbjk\ 2 S{uJ — r2f,fc)/ b ± (cj) is the power spectrum of 
the 6th reservoir, = P J^duj' /tt)^{uj')/ {uj' - 

uj), and xl — X, XR = 0- Here, P means the principal 
value, f b + (u) = 1/(1 + e -ft(1"-Mb)) i s the fermi distri- 
bution function with pbip^b-, and f^(^) = 1 — /b~( w )- 
The control parameters a in this model can be both the 
system parameters, the levels {e^} of the dots and the 
transfer v between the dots, and the reservoir parameters 
{Am Mb}- 

In this model within the RWA, we analytically obtain 
the BSN vector potential (see Appendix IB"! for the deriva- 
tion) : 

(m dp ss \ _\p r Lj (ujj) d f J2b r bj {Wj ) ft (oJj ) \ 

(17) 

where Tjfa) = T,b T bj(^j), and r bj (uj) = + 
^bj&j) 1S t ne spectral function of the 6th reservoir. Also 
we can calculate the BSN curvature as 

TP 

2-f da m \ Tjfa) J 8a n \ Tjiuj) 1 ' 

(18) 

From Eq. (|18[) , we find that no net excess number of 
electrons flow per cycle if only the reservoir parameters 
(/?Lj Ml, /3r, pr) are modulated with the system parame- 
ters fixed. This is because the spectral function is written 
as T b j(uij) = 2-7T J2k Wbjk\ 2 8(uj - tlbk), so that T^j/Tj 
in Eq. (fT8|) is independent of (/3l, Pl, Pr, Pr)- We note 
that this result is characteristic of fermion systems. In- 
deed, in a model of single two-level system connected 
to two bosonic heat reservoirs, there exists heat pump 
by cyclic modulations of the temperatures of the two 
reservoirs 2 ^. This difference between the results for the 
fermion and boson reservoirs comes from the particle 
statistics, which leads that (pc + cc*)b is independent 
of (dependent on) the reservoir parameters for fermion 
(boson) reservoir. Because and 3>j^. are respectively 

proportional to (c'c)b and (cc')&, the particle statistics 
determines whether Tbj = $i depends on the reser- 
voir parameters. 

Before closing this subsection, we make a remark on 
the validity of the RWA on the transport systems. The 
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GQME gives the same result either with or without the 
RWA, as far as the transport between the system and the 
reservoirs is studied, whereas it is known that the internal 
current in the system vanishes in nonequilibrium steady 
states under the RWA— 9 . In Appendix [C] we analytically 
show that the unit-time cumulant generating function 
Aq (a) of the quantity transferred from the reservoirs to 
the system in the steady state for fixed a is equivalent 
between the GQMEs within and without the RWA. We 
also confirm numerically that the results of the adiabatic 
pump within the RWA quantitatively agree with those 
without the RWA in the next subsection. 



B. Non-interacting electron in double quantum dot 
without RWA 

Next we consider a non-interacting double quantum 
dot system coupled to two reservoirs, as illustrated in 
Fig. [2] (a). We assume the wide band limit: T b j(uj) = 
r^. = const. In this subsection we use three different 
methods for calculating (Z\g)° x in modulations of the con- 
trol parameters of the model. In the first method (de- 
noted by RWA), we apply our results within the RWA 
[Eqs. (fTT)) and (|18l) ]. In the second one (denoted by Non- 
RWAl), we numerically solve the eigenvalue problem of 
the GQME generator K x without RWA, and use Eq. ([12]). 
In the third method (denoted by NonRWA2), we numer- 
ically solve the time evolution differential equation of the 
GQME without RWA, and use S T (x) = Tr s ^(r). 

In Fig.[2](b), we plot the excess electron number (Aq)° x 
transferred from the reservoir L to the system for non- 
cyclic modulations of and /zr along the curves 
illustrated in the inset of Fig. [5] (b). We plot the depen- 
dence of (Aq)^ on the final value of the right reservoir 
chemical potential /^ R n of the modulations. We see that 
all the results agree within the numerical precision. This 
implies that the no-pumping condition described below 
Eq. ([18]) within the RWA still holds without the RWA. 
We have also confirmed that the absolute value of the 
BSN curvature computed without the RWA is less than 
10~ 6 in the space of the reservoir parameters, which is 
zero within the numerical precision. 

In Fig. [2] (c), we plot the BSN curvature F £1£2 calcu- 
lated by the method of NonRWAl as a function of E\ and 
£2- We see that the curvature takes the non-zero values 
in this case. Therefore we have a finite geometrical pump 
for the slow periodic modulation of the dot levels (system 
parameters). We note that this result of the BSN curva- 
ture also agrees with the RWA result given by Eq. (TT81) . 
although not shown in the figure. 

We also calculate (Aq)° K for the cyclic process depicted 
in Fig. [1(c). In Fig. [1(d), we plot (A?)° x calculated by 
NonRWA2 as a function of the cycle period r c for various 
values of the amplitude r° ■ of the spectral function. We 
see that for large r c the asymptotic results by NonRWA2 
agree with the results by RWA and NonRWAl. This sup- 
ports the validity of the adiabatic approximation used in 
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FIG. 2. (Color online) (a) Quantum double dot coupled to 
reservoirs, (b) Excess electron transfer {Aq)^ when the reser- 
voir parameters /9l, Ml, and ^r are slowly varied along the 
curves (A, B, and C) shown in the inset, where /3r and the sys- 
tem parameters are fixed (/3r = 1 [1/K], £i = £2 = 0.5 [meV], 
and v — 0.2 [meV]). The initial values of the parameters 
are /3 L ni = 1 [1/K] and ^ = ^ = 0.2 [meV]. The de- 
pendence on the final value fi^ n is plotted, where /3^ n — /3 L nl 
and /Xl" = Ml"- The calculations are performed by RWA and 
NonRWAl. (c) BSN curvature F ei e 2 as a function of £1 and 
£2, calculated by NonRWAl. The reservoir parameters are 
set to /3 L = /3 R = 1 [1/K] and ^ L = MR = 0.2 [meV]. (d) 
Excess electron transfer (Aq)^ for a cyclic process along the 
circle depicted in (c). The horizontal axis is the period r c of 
the cyclic process. The calculations are performed by RWA, 
NonRWAl, and NonRWA2. The results by NonRWA2 are 
plotted for various values of the amplitude of the reser- 
voir spectral function, while it is fixed to 0.001 [meV] for the 
other methods. 



deriving Eq. (|12p for slow modulations. The character- 
istic time for the validity of the adiabatic approxima- 
tion becomes shorter as 1^ increases, which implies that 
t c 3> h/Tjjj is an adiabatic condition. 

We note that all the results by RWA and NonRWAl 
agree with each other not only qualitatively but also 
quantitatively. This implies that the rotating wave ap- 
proximation is valid in discussing transport between the 
system and the reservoirs under slow modulations of the 
parameters. 



C. Interacting electron model 

We next consider an interacting spinless electron sys- 
tem in a double quantum dot. We take the system Hamil- 
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(hiLi2 + U,huji + U). This implies that the positive and 
negative peaks merge for non- interacting system (U = 0), 
and thus the curvature becomes zero in (^l,Mr)~ 

space to achieve the no-pumping condition mentioned be- 
low Eq. (fTS)) . This result also implies that, for U > 0, 
an adiabatic pump can occur even if only the reservoir 
parameters are modulated. 

Indeed, Fig. [3] (e) shows the [/-dependence of the ex- 
cess electron transfer (Aq)° x for cyclic processes [cycle 1 
and 2 depicted in Fig. [3J (c)], where (Aq)°* is non-zero 
for U > 0. Note that the direction of the cycle 2 is op- 
posite to that of the cycle 1. We also note that (Aq)°* 
for the cycle 1 becomes a constant for U > l[meV], be- 
cause the positive peak positions of F MLAtR do not move 
as U increases. In contrast, (Aq) e * for the cycle 2 has 
a peak at U ~ 0.5 [me V] and becomes nearly zero for 
U >1 [meV] . This results from the fact that the negative 
peak positions of F^ L ^ R move away from the cycle 2 as 
U increases. 



FIG. 3. (Color online) (a)-(c) BSN curvature F MlMr for the 
interacting electron system in the double dot, as a function of 
/xl and /xr. The interaction strength is U = 0.25[meV] for (a), 
U = 0.5[meV] for (b), and U = 1.0[meV] for (c). The other 
parameters are /3l = Pn = 1 [1/K], ei = £2 = 0.5 [meV], 
v = 0.2 [meV], and = 0.001 [meV]. (d) Positions of the 
peaks of the BSN curvature, as a function of U. The position 
(/Xl° s \ Mr 551 ) °f the positive peak with /J[° sl > Mr° s1 and that 
(/x£ ega , /x R Gga ) of the negative peaks with /x£ cga > ^ R cga are 
shown. The dashed lines are the mode energies, Hu)i and hui2 , 
of the system Hamiltonian (fiwi < tiw^)- The solid lines are 
huj\ + U and fiu2 + U. (e) Excess electron transfer (Aq)^ 
for cyclic processes along the circles (1 and 2) depicted in 
(c). The horizontal axis is the interaction strength U. The 
direction of the cycle 2 is opposite to that of the cycle 1. 



tonian as 

H s = £ i^i + v (d\d2 + 4 rf "i) + Ud\didld 2 - (19) 

i=l,2 

In this model, an electron in one dot interacts with an 
electron in the other dot. We investigate the excess elec- 
tron transfer (Aq}^ x with the modulation of the chemical 
potentials (/il and /xr). We also assume the wide band 
limit: Tf,j(ui) = = const. 

In Fig. [3] (a)-(c), we plot the BSN curvature i^ LMR 
for various values of the interaction strength U obtained 
from the numerical diagonalization method without the 
RWA. We observe the two positive and two negative 
peaks. In Fig. |3J (d), we plot the position of one of the 
positive peaks and that of one of the negative peaks. 
We see that the positions of the negative peaks move 
as U increases, whereas those of the positive ones do 
not. Moreover we find that these peak positions are lo- 
cated around at the energies necessary to add one elec- 
tron; the approximate position of the positive peaks are 
(ml, Mr) = (fiwij fu^2} and (liW2,Swi), and those of the 
negative ones are (/ji^Mr) = (^1 + U 1 hi02 + U) and 



IV. DISCUSSION AND CONCLUSION 

By using a QME approach, we have derived a geomet- 
rical expressions of the cumulant generating function and 
the average of the pumped (excess) quantity transferred 
from a reservoir to the system under slow modulation of 
control parameters, where the BSN phases and the BSN 
curvature of the QME play crucial roles. 

For non-interacting electrons in quantum dot systems, 
there is no pump current when only the temperatures 
and chemical potentials of the reservoirs are modulated. 
In contrast, for an interacting system, the pump current 
can be observed even in this situation. We note that 
the modulations of only the chemical potentials of the 
reservoirs are required for the pump. This has an advan- 
tage for the control of the pump in experiments, since 
the modulation of chemical potential is easier than that 
of the temperatures. As shown in Fig. [3J (e), for a cyclic 
modulation of the chemical potentials, the pump current 
depends not only on the difference of the chemical po- 
tentials but also on their average. This implies that the 
average number of electrons in the system S is impor- 
tant for the pump. For example, when we modulate the 
chemical potentials as ml(^) = Mc + Mr cos(27ri/T c ) and 
MrM = Mc + Mr sin(27ri/r c ), /i c affects the quantity of 
the pump. This fact may be applicable for switching the 
pump by the change of fi c or the electron density, which 
can be controlled by a gate voltage. 

Since we applied the method of the full counting statis- 
tics, we can calculate the fluctuation of the pump vari- 
able. It is a future issue to analyse the detailed properties 
of the fluctuation in the adiabatic pump. It is also in- 
teresting to investigate the relation between the present 
geometrical expression of the adiabatic pump based on 
the QME and the conventional geometrical expressions 
based on the scattering theory^—. The investigations of 
non-adiabatic pump, non-Markovian situation, and spin 
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effect are also future issues. 
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Appendix A: Generalized quantum master equation 
without and within the RWA 

We here derive the concrete form of the GQME (JS])-© 
without and within the RWA. 

We start from the total Hamiltonian (system plus 
reservoirs) given as 

fftot - H s + Hi, + u J2 H Sb - (Al) 

b b 

Here, for simplicity, we assume that the reservoir Hamil- 
tonian H b and the coupling Hamiltonian Hsb between 
the system and the reservoir are respectively written as 

Hb = Y Kl bk c\ k c b k, (A2) 
fc 

H Sb = J2( Vbk& l £bk + V bk< d lk)> (A3) 
fe 

where is a single-particle mode annihilation operator 
in the system S, ib is the index of the system mode that 



couples to the 6th reservoir, and c b k is the fcth mode 
annihilation operator in the 6th reservoir. We denote the 
eigenenergy of the system Hamiltonian Hs as E x , and the 
corresponding energy eigenstate as \E X ). We also assume 
that all the eigenenergies of Hg are non-degenerate. 

We consider the quantity Q = J2bJ2k 1bkc\ k c b k to be 
counted, and define the current of Q from the reservoirs 
to the system S as positive. 

For the derivation of the GQME, it is convenient to 
introduce the eigenoperatorsS? from : 

«h S) = \ E * - huj s){E x - tkJ s \a h \E x )(E x \, (A4) 
"1 (WS) = E \ E * + h "s)(E x + hw s \a\ b \E x ){E x \. (A5) 

Ex 

Then the modified coupling Hamiltonian in the interac- 
tion picture is written as 

Hg b (t) = e - i xQ/2 e -(Hs+H b )t/ihjjX be (H s +H,,)t/ih e ixQ/2 
k us 

+ Vjjo^^e-*™*/^-^-"")*) . (A6) 

We assume that the initial state of the total system is 
written as /O to t(0) = Po^Pies-, where po is a initial state of 
the system S, p ICS = (g) b pf, and pf = e'^ 6b ~^"^ /Z b 
is the grand-canonical state of the 6th reservoir. Then 
substituting Eq. (|A6|) into Eq. ©, we obtain the GQME 



b OJgtJg 

+ $ b -K){«r s) a^ ) /5 x (0 + pHt)al^a^ - e~™^ + a^ k) p*(t)a£*>) } 

+ if t(u' s ){-a^4^px(t) + pHt)a^al^} + .^^{-^W^ + ^V(*)< (Ws) } 

(A7) 



where we used 



Here, q(fl b k) = qbk, and 

$f («) = Y, toWbk - uWbkfffc, (A9) 



dt'e iut ' 



(A8) 



, i , , , .) D \^ bk \ 2 fbk. 

Q 6fc - U3 ' 



(A10) 
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* b ± (^; X )=2^P 



\ v bk\ 2 fbk c ±t X g(n bk ) ^ {Ml) 

&bk — W 



fbk = ^b{pbcl k c bk } = - 
fbk = TT b{pbC bk cl k } = I - f, 



1 



(A12) 



e l3 b (hn bk -fj, b ) ' 

(A13) 



Equation (IA7|1 is the concrete form of the GQME without 
the RWA. 

When we transform Eq. (|A7j) in the interaction pic- 
ture, we see that rapidly oscillating terms proportional 
to exp[±i(ws— Wg)t] appear. In the RWA we neglect these 
terms2£. Thus we obtain the GQME with the RWA by 
leaving only the terms with u' s = ui$ in Eq. (|A7[) . 

For the non-interacting models in Sees. IIIII A and 
B, because the eigenoperators are the mode operators 
themselves, the GQMEs for these models, in particular 
Eq. (fTB) for the RWA case, are derived. 



Appendix B: Derivation of Eq. (|17[) 

In the model in Sec. IIIII A within the RWA, the eigen- 
values and the eigenvectors of /C x can be decomposed into 
those of fC xJ . That is, A£ = £j A^, £g = <g) j 1^, and 
/5q = . where Aq ^ is the eigenvalue of IC x j with 

maximum real part, and £qj and /5g ■ are respectively 
the corresponding left and right eigenvectors, which are 
operators on the jth mode Hilbert space. 

When we represent the left and right eigenvectors in 
the basis of the number states (denoted by \0j) and | lj) ) 
of a,j such that aj\0j) = and \lj) = aAOj), we can 
show that (mjl^o |(1 — m)j) — (m,j\p^ -|(1 — m)j) = 



(to = 0, 1), and 



i 

«i(x) 



(Bl) 
(B2) 



where 



E b K-K)-^-K)) + /P7 

2(*y(u i )e-*x + $^( Wi )) 
D i =r i 2 -4(l-e^)$+(a; i )*R,(w 3 -) 



u i(x) 
w 3 (x) 



-4(l- e -^)$- j K)<I>+ J (^), 



r j = Eb T bj, and r hj = + ^ bj {u>j). From the 

normalization condition for x — 0, TrsPp^- = 1, we have 
Cj(0) = £ h $y(wj)/rj. Thus we obtain the BSN vector 
potential: 



Trc £ 







da. 



d(ix) 



a(c 3 -(0H-(0)) 
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(B3) 



which becomes the desired result after straightforward 
calculation. 



Appendix C: Equivalence of the unit-time cumulant 
generating functions with and without the RWA 

1. Matrix representation of GQME generator 

To show the equivalence, we introduce the matrix rep- 
resentation of JC X by using the eigenstates \E X ) of Hs- 
the (y'y,x'x) matrix element is given as (lC x )y' v ,x'x 
Tr 



where xx' (yy') is the 



(\E y ,)(E y \) f (K x \E x ,)(E x \) 

index for the column (row) of the matrix. 

In this representation, we can show that within the 
= (X^ WA )„ Wl x'« i: •<•' ' x 
is a block di- 



This implies that K^ WA 



RWA (lC™ A ) y ^ xx 

and y' ^ y. This „. n , N 

agonal matrix that is composed of {(fc^^yy^x} and 
{(/C* WA W^} with + x and y' ^ y. We also 
note that a relation holds between the matrices of the 
generators without and with the RWA: (/C^ WA )j,j, ia;K — 
(fc x )yy,xx- 



2. Equivalence of Aq with and without RWA 

As is mentioned in Sec. [Til B, the unit-time cumulant 
generating function in a steady state is given by the eigen- 
value Aq of the generator IC X with maximum real part. 

Within the RWA, Aq is determined by the eigenvalue 
of {{1C* WA ) yy , xx }, one of the blocks of {{K.™ A ) W , XX ,}. 

Without the RWA, the eigenvalues of K, x is determined 
by a perturbation theory with respect to v = u 2 . From 
Eq. (|AT[> . we see that the unperturbed part of K x is 
—{{E x —E x i)/ih}8E y ,E m $E i,e x ,, and is diagonal. There- 
fore the unperturbed eigenvalue is —(E x — E x >)/ih. This 
implies that the eigenvalue of zero has d-fold degener- 
acy, where d = dimHs- Thus by the perturbation the- 
ory for degenerate case, the first order eigenvalue is de- 
termined by the eigenvalue equation for the matrix in 
the degenerate subspace, i.e., {{lC x )yy,xx}- Furthermore, 
since (fC x )y y ,xx = {lC s ^ WA ) yy _ xx , the first order eigenval- 
ues are equivalent to those of {tC™ ) yy , X x- Therefore Ag 
without the RWA is equivalent to Ag within the RWA 
in 0(v). This verification of the equivalence is sufficient 
since the master equation is valid up to 0{v). 
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